Tutorial on how to build non-Markovian dynamic models from molecular dynamics simulations for studying protein conformational changes.

Protein conformational changes play crucial roles in their biological functions. In recent years, the Markov State Model (MSM) constructed from extensive Molecular Dynamics (MD) simulations has emerged as a powerful tool for modeling complex protein conformational changes. In MSMs, dynamics are modeled as a sequence of Markovian transitions among metastable conformational states at discrete time intervals (called lag time). A major challenge for MSMs is that the lag time must be long enough to allow transitions among states to become memoryless (or Markovian). However, this lag time is constrained by the length of individual MD simulations available to track these transitions. To address this challenge, we have recently developed Generalized Master Equation (GME)-based approaches, encoding non-Markovian dynamics using a time-dependent memory kernel. In this Tutorial, we introduce the theory behind two recently developed GME-based non-Markovian dynamic models: the quasi-Markov State Model (qMSM) and the Integrative Generalized Master Equation (IGME). We subsequently outline the procedures for constructing these models and provide a step-by-step tutorial on applying qMSM and IGME to study two peptide systems: alanine dipeptide and villin headpiece. This Tutorial is available at https://github.com/xuhuihuang/GME_tutorials. The protocols detailed in this Tutorial aim to be accessible for non-experts interested in studying the biomolecular dynamics using these non-Markovian dynamic models.

Submillisecond Atomistic Molecular Dynamics Simulations Reveal Hydrogen Bond-Driven Diffusion of a Guest Peptide in Protein-RNA Condensate.

Liquid-liquid phase separation mediated by proteins and/or nucleic acids is believed to underlie the formation of many distinct condensed phases, or membraneless organelles, within living cells. These condensates have been proposed to orchestrate a variety of important processes. Despite recent advances, the interactions that regulate the dynamics of molecules within a condensate remain poorly understood. We performed accumulated 564.7 µs all-atom molecular dynamics (MD) simulations (system size ∼200k atoms) of model condensates formed by a scaffold RNA oligomer and a scaffold peptide rich in arginine (Arg). These model condensates contained one of three possible guest peptides: the scaffold peptide itself or a variant in which six Arg residues were replaced by lysine (Lys) or asymmetric dimethyl arginine (ADMA). We found that the Arg-rich peptide can form the largest number of hydrogen bonds and bind the strongest to the scaffold RNA in the condensate, relative to the Lys- and ADMA-rich peptides. Our MD simulations also showed that the Arg-rich peptide diffused more slowly in the condensate relative to the other two guest peptides, which is consistent with a recent fluorescence microscopy study. There was no significant increase in the number of cation-π interactions between the Arg-rich peptide and the scaffold RNA compared to the Lys-rich and ADMA-rich peptides. Our results indicate that hydrogen bonds between the peptides and the RNA backbone, rather than cation-π interactions, play a major role in regulating peptide diffusion in the condensate.

Integrative generalized master equation: A method to study long-timescale biomolecular dynamics via the integrals of memory kernels.

The generalized master equation (GME) provides a powerful approach to study biomolecular dynamics via non-Markovian dynamic models built from molecular dynamics (MD) simulations. Previously, we have implemented the GME, namely the quasi Markov State Model (qMSM), where we explicitly calculate the memory kernel and propagate dynamics using a discretized GME. qMSM can be constructed with much shorter MD trajectories than the MSM. However, since qMSM needs to explicitly compute the time-dependent memory kernels, it is heavily affected by the numerical fluctuations of simulation data when applied to study biomolecular conformational changes. This can lead to numerical instability of predicted long-time dynamics, greatly limiting the applicability of qMSM in complicated biomolecules. We present a new method, the Integrative GME (IGME), in which we analytically solve the GME under the condition when the memory kernels have decayed to zero. Our IGME overcomes the challenges of the qMSM by using the time integrations of memory kernels, thereby avoiding the numerical instability caused by explicit computation of time-dependent memory kernels. Using our solutions of the GME, we have developed a new approach to compute long-time dynamics based on MD simulations in a numerically stable, accurate and efficient way. To demonstrate its effectiveness, we have applied the IGME in three biomolecules: the alanine dipeptide, FIP35 WW-domain, and Taq RNA polymerase. In each system, the IGME achieves significantly smaller fluctuations for both memory kernels and long-time dynamics compared to the qMSM. We anticipate that the IGME can be widely applied to investigate biomolecular conformational changes.

Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently.

Conformational changes underpin function and encode complex biomolecular mechanisms. Gaining atomic-level detail of how such changes occur has the potential to reveal these mechanisms and is of critical importance in identifying drug targets, facilitating rational drug design, and enabling bioengineering applications. While the past two decades have brought Markov state model techniques to the point where practitioners can regularly use them to glimpse the long-time dynamics of slow conformations in complex systems, many systems are still beyond their reach. In this Perspective, we discuss how including memory (i.e., non-Markovian effects) can reduce the computational cost to predict the long-time dynamics in these complex systems by orders of magnitude and with greater accuracy and resolution than state-of-the-art Markov state models. We illustrate how memory lies at the heart of successful and promising techniques, ranging from the Fokker-Planck and generalized Langevin equations to deep-learning recurrent neural networks and generalized master equations. We delineate how these techniques work, identify insights that they can offer in biomolecular systems, and discuss their advantages and disadvantages in practical settings. We show how generalized master equations can enable the investigation of, for example, the gate-opening process in RNA polymerase II and demonstrate how our recent advances tame the deleterious influence of statistical underconvergence of the molecular dynamics simulations used to parameterize these techniques. This represents a significant leap forward that will enable our memory-based techniques to interrogate systems that are currently beyond the reach of even the best Markov state models. We conclude by discussing some current challenges and future prospects for how exploiting memory will open the door to many exciting opportunities.

EPISOL: A software package with expanded functions to perform 3D-RISM calculations for the solvation of chemical and biological molecules.

Integral equation theory (IET) provides an effective solvation model for chemical and biological systems that balances computational efficiency and accuracy. We present a new software package, the expanded package for IET-based solvation (EPISOL), that performs 3D-reference interaction site model (3D-RISM) calculations to obtain the solvation structure and free energies of solute molecules in different solvents. In EPISOL, we have implemented 22 different closures, multiple free energy functionals, and new variations of 3D-RISM theory, including the recent hydrophobicity-induced density inhomogeneity (HI) theory for hydrophobic solutes and ion-dipole correction (IDC) theory for negatively charged solutes. To speed up the convergence and enhance the stability of the self-consistent iterations, we have introduced several numerical schemes in EPISOL, including a newly developed dynamic mixing approach. We show that these schemes have significantly reduced the failure rate of 3D-RISM calculations compared to AMBER-RISM software. EPISOL consists of both a user-friendly graphic interface and a kernel library that allows users to call its routines and adapt them to other programs. EPISOL is compatible with the force-field and coordinate files from both AMBER and GROMACS simulation packages. Moreover, EPISOL is equipped with an internal memory control to efficiently manage the use of physical memory, making it suitable for performing calculations on large biomolecules. We demonstrate that EPISOL can efficiently and accurately calculate solvation density distributions around various solute molecules (including a protein chaperone consisting of 120,715 atoms) and obtain solvent free energy for a wide range of organic compounds. We expect that EPISOL can be widely applied as a solvation model for chemical and biological systems. EPISOL is available at https://github.com/EPISOLrelease/EPISOL.

Building insightful, memory-enriched models to capture long-time biochemical processes from short-time simulations.

The ability to predict and understand complex molecular motions occurring over diverse timescales ranging from picoseconds to seconds and even hours in biological systems remains one of the largest challenges to chemical theory. Markov state models (MSMs), which provide a memoryless description of the transitions between different states of a biochemical system, have provided numerous important physically transparent insights into biological function. However, constructing these models often necessitates performing extremely long molecular simulations to converge the rates. Here, we show that by incorporating memory via the time-convolutionless generalized master equation (TCL-GME) one can build a theoretically transparent and physically intuitive memory-enriched model of biochemical processes with up to a three order of magnitude reduction in the simulation data required while also providing a higher temporal resolution. We derive the conditions under which the TCL-GME provides a more efficient means to capture slow dynamics than MSMs and rigorously prove when the two provide equally valid and efficient descriptions of the slow configurational dynamics. We further introduce a simple averaging procedure that enables our TCL-GME approach to quickly converge and accurately predict long-time dynamics even when parameterized with noisy reference data arising from short trajectories. We illustrate the advantages of the TCL-GME using alanine dipeptide, the human argonaute complex, and FiP35 WW domain.

The Ion-Dipole Correction of the 3DRISM Solvation Model to Accurately Compute Water Distributions around Negatively Charged Biomolecules.

The 3D reference interaction site model (3DRISM) provides an efficient grid-based solvation model to compute the structural and thermodynamic properties of biomolecules in aqueous solutions. However, it remains challenging for existing 3DRISM methods to correctly predict water distributions around negatively charged solute molecules. In this paper, we first show that this challenge is mainly due to the orientation of water molecules in the first solvation shell of the negatively charged solute molecules. To properly consider this orientational preference, position-dependent two-body intramolecular correlations of solvent need to be included in the 3DRISM theory, but direct evaluations of these position-dependent two-body intramolecular correlations remain numerically intractable. To address this challenge, we introduce the Ion-Dipole Correction (IDC) to the 3DRISM theory, in which we incorporate the orientation preference of water molecules via an additional solute-solvent interaction term (i.e., the ion-dipole interaction) while keeping the formulism of the 3DRISM equation unchanged. We prove that this newly introduced IDC term is equivalent to an effective direct correlation function which can effectively consider the orientation effect that arises from position dependent two-body correlations. We first quantitatively validate our 3DRISM-IDC theory combined with the PSE3 closure on Cl-, [ClO]- (a two-site anion), and [NO2]- (a three-site anion). For all three anions, we show that our 3DRISM-IDC theory significantly outperforms the 3DRISM theory in accurately predicting the solvation structures in comparison to MD simulations, including RDFs and 3D water distributions. Furthermore, we have also demonstrated that the 3DRISM-IDC can improve the accuracy of hydration free-energy calculation for Cl-. We further demonstrate that our 3DRISM-IDC theory yields significant improvements over the 3DRISM theory when applied to compute the solvation structures for various negatively charged solute molecules, including adenosine triphosphate (ATP), a short peptide containing 19 residues, a DNA hairpin containing 24 nucleotides, and a riboswitch RNA molecule with 77 nucleotides. We expect that our 3DRISM-IDC-PSE3 solvation model holds great promise to be widely applied to study solvation properties for nucleic acids and other biomolecules containing negatively charged functional groups.

Elucidation of the key role of Pt···Pt interactions in the directional self-assembly of platinum(II) complexes.

Here, we report the use of an amphiphilic Pt(II) complex, K[Pt{(O3SCH2CH2CH2)2bzimpy}Cl] (PtB), as a model to elucidate the key role of Pt···Pt interactions in directing self-assembly by combining temperature-dependent ultraviolet-visible (UV-Vis) spectroscopy, stopped-flow kinetic experiments, quantum mechanics (QM) calculations, and molecular dynamics (MD) simulations. Interestingly, we found that the self-assembly mechanism of PtB in aqueous solution follows a nucleation-free isodesmic model, as revealed by the temperature-dependent UV-Vis experiments. In contrast, a cooperative growth is found for the self-assembly of PtB in acetone­water (7:1, vol/vol) solution, which is further verified by the stopped-flow experiments, which clearly indicates the existence of a nucleation phase in the acetone­water (7:1, vol/vol) solution. To reveal the underlying reasons and driving forces for these self-assembly processes, we performed QM calculations and show that the Pt···Pt interactions arising from the interaction between the pz and dz2 orbitals play a crucial role in determining the formation of ordered self-assembled structures. In subsequent oligomer MD simulations, we demonstrate that this directional Pt···Pt interaction can indeed facilitate the formation of linear structures packed in a helix-like fashion. Our results suggest that the self-assembly of PtB in acetone­water (7:1, vol/vol) solution is predominantly driven by the directional noncovalent Pt···Pt interaction, leading to the cooperative growth and the formation of fibrous nanostructures. On the contrary, the self-assembly in aqueous solution forms spherical nanostructures of PtB, which is primarily due to the predominant contribution from the less directional hydrophobic interactions over the directional Pt···Pt and π−π interactions that result in an isodesmic growth.

RPnet: a reverse-projection-based neural network for coarse-graining metastable conformational states for protein dynamics.

The Markov State Model (MSM) is a powerful tool for modeling long timescale dynamics based on numerous short molecular dynamics (MD) simulation trajectories, which makes it a useful tool for elucidating the conformational changes of biological macromolecules. By partitioning the phase space into discretized states and estimating the probabilities of inter-state transitions based on short MD trajectories, one can construct a kinetic network model that could be used to extrapolate long-timescale kinetics if the Markovian condition is met. However, meeting the Markovian condition often requires hundreds or even thousands of states (microstates), which greatly hinders the comprehension of the conformational dynamics of complex biomolecules. Kinetic lumping algorithms can coarse grain numerous microstates into a handful of metastable states (macrostates), which would greatly facilitate the elucidation of biological mechanisms. In this work, we have developed a reverse-projection-based neural network (RPnet) to lump microstates into macrostates, by making use of a physics-based loss function that is based on the projection operator framework of conformational dynamics. By recognizing that microstate and macrostate transition modes can be related through a projection process, we have developed a reverse-projection scheme to directly compare the microstate and macrostate dynamics. Based on this reverse-projection scheme, we designed a loss function that allows the effective assessment of the quality of a given kinetic lumping. We then make use of a neural network to efficiently minimize this loss function to obtain an optimized set of macrostates. We have demonstrated the power of our RPnet in analyzing the dynamics of a numerical 2D potential, alanine dipeptide, and the clamp opening of an RNA polymerase. In all these systems, we have illustrated that our method could yield comparable or better results than competing methods in terms of state partitioning and reproduction of slow dynamics. We expect that our RPnet holds promise in analyzing the conformational dynamics of biological macromolecules.

Critical role of backbone coordination in the mRNA recognition by RNA induced silencing complex.

Despite its functional importance, the molecular mechanism underlying target mRNA recognition by Argonaute (Ago) remains largely elusive. Based on extensive all-atom molecular dynamics simulations, we constructed quasi-Markov State Model (qMSM) to reveal the dynamics during recognition at position 6-7 in the seed region of human Argonaute 2 (hAgo2). Interestingly, we found that the slowest mode of motion therein is not the gRNA-target base-pairing, but the coordination of the target phosphate groups with a set of positively charged residues of hAgo2. Moreover, the ability of Helix-7 to approach the PIWI and MID domains was found to reduce the effective volume accessible to the target mRNA and therefore facilitate both the backbone coordination and base-pair formation. Further mutant simulations revealed that alanine mutation of the D358 residue on Helix-7 enhanced a trap state to slow down the loading of target mRNA. Similar trap state was also observed when wobble pairs were introduced in g6 and g7, indicating the role of Helix-7 in suppressing non-canonical base-paring. Our study pointed to a general mechanism for mRNA recognition by eukaryotic Agos and demonstrated the promise of qMSM in investigating complex conformational changes of biomolecular systems.

Markov State Models to Study the Functional Dynamics of Proteins in the Wake of Machine Learning.

Markov state models (MSMs) based on molecular dynamics (MD) simulations are routinely employed to study protein folding, however, their application to functional conformational changes of biomolecules is still limited. In the past few years, the field of computational chemistry has experienced a surge of advancements stemming from machine learning algorithms, and MSMs have not been left out. Unlike global processes, such as protein folding, the application of MSMs to functional conformational changes is challenging because they mostly consist of localized structural transitions. Therefore, it is critical to properly select a subset of structural features that can describe the slowest dynamics of these functional conformational changes. To address this challenge, we recommend several automatic feature selection methods such as Spectral-OASIS. To identify states in MSMs, the chosen features can be subject to dimensionality reduction methods such as TICA or deep learning based VAMPNets to project MD conformations onto a few collective variables for subsequent clustering. Another challenge for the application of MSMs to the study of functional conformational changes is the ability to comprehend their biophysical mechanisms, as MSMs built for these processes often require a large number of states. We recommend the recently developed quasi-MSMs (qMSMs) to address this issue. Compared to MSMs, qMSMs encode the non-Markovian dynamics via the generalized master equation and can significantly reduce the number of states. As a result, qMSMs can be built with a handful of states to facilitate the interpretation of functional conformational changes. In the wake of machine learning, we believe that the rapid advancement in the MSM methodology will lead to their wider application in studying functional conformational changes of biomolecules.

Entropy of stapled peptide inhibitors in free state is the major contributor to the improvement of binding affinity with the GK domain.

Stapled peptides are promising protein-protein interaction (PPI) inhibitors that can increase the binding potency. Different from small-molecule inhibitors in which the binding mainly depends on energetic interactions with their protein targets, the binding of stapled peptides has long been suggested to be benefited from entropy. However, it remains challenging to reveal the molecular features that lead to this entropy gain, which could originate from the stabilization of the stapled peptide in solution or from the increased flexibility of the complex upon binding. This hinders the rational design of stapled peptides as PPI inhibitors. Using the guanylate kinase (GK) domain of the postsynaptic density protein 95 (PSD-95) as the target, we quantified the enthalpic and entropic contributions by combining isothermal titration calorimetry (ITC), X-ray crystallography, and free energy calculations based on all-atom molecular dynamics (MD) simulations. We successfully designed a stapled peptide inhibitor (staple 1) of the PSD-95 GK domain that led to a 25-fold increase in the binding affinity (from tens of µMs to 1.36 µM) with high cell permeability. We showed that entropy indeed greatly enhanced the binding affinity and the entropy gain was mainly due to the constrained-helix structure of the stapled peptide in solution (free state). Based on staple 1, we further designed two other stapled peptides (staple 2 and 3), which exerted even larger entropy gains compared to staple 1 because of their more flexible bound complexes (bound state). However, for staple 2 and 3, the overall binding affinities were not improved, as the loose binding in their bound states led to an enthalpic loss that largely compensated the excess entropy gain. Our work suggests that increasing the stability of the stapled peptide in free solution is an effective strategy for the rational design of stapled peptides as PPI inhibitors.

Role of bacterial RNA polymerase gate opening dynamics in DNA loading and antibiotics inhibition elucidated by quasi-Markov State Model.

To initiate transcription, the holoenzyme (RNA polymerase [RNAP] in complex with σ factor) loads the promoter DNA via the flexible loading gate created by the clamp and ß-lobe, yet their roles in DNA loading have not been characterized. We used a quasi-Markov State Model (qMSM) built from extensive molecular dynamics simulations to elucidate the dynamics of Thermus aquaticus holoenzyme's gate opening. We showed that during gate opening, ß-lobe oscillates four orders of magnitude faster than the clamp, whose opening depends on the Switch 2's structure. Myxopyronin, an antibiotic that binds to Switch 2, was shown to undergo a conformational selection mechanism to inhibit clamp opening. Importantly, we reveal a critical but undiscovered role of ß-lobe, whose opening is sufficient for DNA loading even when the clamp is partially closed. These findings open the opportunity for the development of antibiotics targeting ß-lobe of RNAP. Finally, we have shown that our qMSMs, which encode non-Markovian dynamics based on the generalized master equation formalism, hold great potential to be widely applied to study biomolecular dynamics.

On the advantages of exploiting memory in Markov state models for biomolecular dynamics.

Biomolecular dynamics play an important role in numerous biological processes. Markov State Models (MSMs) provide a powerful approach to study these dynamic processes by predicting long time scale dynamics based on many short molecular dynamics (MD) simulations. In an MSM, protein dynamics are modeled as a kinetic process consisting of a series of Markovian transitions between different conformational states at discrete time intervals (called "lag time"). To achieve this, a master equation must be constructed with a sufficiently long lag time to allow interstate transitions to become truly Markovian. This imposes a major challenge for MSM studies of proteins since the lag time is bound by the length of relatively short MD simulations available to estimate the frequency of transitions. Here, we show how one can employ the generalized master equation formalism to obtain an exact description of protein conformational dynamics both at short and long time scales without the time resolution restrictions imposed by the MSM lag time. Using a simple kinetic model, alanine dipeptide, and WW domain, we demonstrate that it is possible to construct these quasi-Markov State Models (qMSMs) using MD simulations that are 5-10 times shorter than those required by MSMs. These qMSMs only contain a handful of metastable states and, thus, can greatly facilitate the interpretation of mechanisms associated with protein dynamics. A qMSM opens the door to the study of conformational changes of complex biomolecules where a Markovian model with a few states is often difficult to construct due to the limited length of available MD simulations.

TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromolecules.

Locating the minimum free energy paths (MFEPs) between two conformational states is among the most important tasks of biomolecular simulations. For example, knowledge of the MFEP is critical for focusing the effort of unbiased simulations that are used for the construction of Markov state models to the biologically relevant regions of the system. Typically, existing path searching methods perform local sampling around the path nodes in a pre-selected collective variable (CV) space to allow a gradual downhill evolution of the path toward the MFEP. Despite the wide application of such a strategy, the gradual path evolution and the non-trivial a priori choice of CVs are also limiting its overall efficiency and automation. Here we demonstrate that non-local perpendicular sampling can be pursued to accelerate the search, provided that all nodes are reordered thereafter via a traveling-salesman scheme. Moreover, path-CVs can be computed on-the-fly and used as a coordinate system, minimizing the necessary prior knowledge about the system. Our traveling-salesman based automated path searching method achieves a 5-8 times speedup over the string method with swarms-of-trajectories for two peptide systems in vacuum and solution, making it a promising method for obtaining initial pathways when investigating functional conformational changes between a pair of structures.

Real-time monitoring of hydrophobic aggregation reveals a critical role of cooperativity in hydrophobic effect.

The hydrophobic interaction drives nonpolar solutes to aggregate in aqueous solution, and hence plays a critical role in many fundamental processes in nature. An important property intrinsic to hydrophobic interaction is its cooperative nature, which is originated from the collective motions of water hydrogen bond networks surrounding hydrophobic solutes. This property is widely believed to enhance the formation of hydrophobic core in proteins. However, cooperativity in hydrophobic interactions has not been successfully characterized by experiments. Here, we quantify cooperativity in hydrophobic interactions by real-time monitoring the aggregation of hydrophobic solute (hexaphenylsilole, HPS) in a microfluidic mixer. We show that association of a HPS molecule to its aggregate in water occurs at sub-microsecond, and the free energy change is -5.8 to -13.6 kcal mol-1. Most strikingly, we discover that cooperativity constitutes up to 40% of this free energy. Our results provide quantitative evidence for the critical role of cooperativity in hydrophobic interactions.

Reference interaction site model with hydrophobicity induced density inhomogeneity: An analytical theory to compute solvation properties of large hydrophobic solutes in the mixture of polyatomic solvent molecules.

Reference interaction site model (RISM) has recently become a popular approach in the study of thermodynamical and structural properties of the solvent around macromolecules. On the other hand, it was widely suggested that there exists water density depletion around large hydrophobic solutes (>1 nm), and this may pose a great challenge to the RISM theory. In this paper, we develop a new analytical theory, the Reference Interaction Site Model with Hydrophobicity induced density Inhomogeneity (RISM-HI), to compute solvent radial distribution function (RDF) around large hydrophobic solute in water as well as its mixture with other polyatomic organic solvents. To achieve this, we have explicitly considered the density inhomogeneity at the solute-solvent interface using the framework of the Yvon-Born-Green hierarchy, and the RISM theory is used to obtain the solute-solvent pair correlation. In order to efficiently solve the relevant equations while maintaining reasonable accuracy, we have also developed a new closure called the D2 closure. With this new theory, the solvent RDFs around a large hydrophobic particle in water and different water-acetonitrile mixtures could be computed, which agree well with the results of the molecular dynamics simulations. Furthermore, we show that our RISM-HI theory can also efficiently compute the solvation free energy of solute with a wide range of hydrophobicity in various water-acetonitrile solvent mixtures with a reasonable accuracy. We anticipate that our theory could be widely applied to compute the thermodynamic and structural properties for the solvation of hydrophobic solute.

Interactions of a water-soluble fullerene derivative with amyloid-ß protofibrils: dynamics, binding mechanism, and the resulting salt-bridge disruption.

Alzheimer's disease (AD) is associated with the pathological self-assembly of amyloid-ß (Aß) peptides into ß-sheet-rich oligomers and insoluble amyloid fibrils. Experimental studies reported that 1,2-(dimethoxymethano)fullerene (DMF), a water-soluble fullerene derivative, inhibits strongly Aß peptide aggregation at the early stage. However, the interaction and binding mechanisms are not well understood. In this study, we have investigated the detailed interaction of a DMF molecule with a fibrillar hexamer of full-length Aß42 and the resulting structural alterations by performing multiple all-atom explicit solvent molecular dynamics (MD) simulations. Starting from different initial states with a minimum distance of 2 nm between the DMF and the Aß protofibril, our MD simulations show that the DMF binds to the Aß protofibril via both slow and fast binding processes. Three dominant binding sites are identified: the central hydrophobic core (CHC) site (17LVFFA21), the turn site (27NKGAI31), and the C-terminal ß-sheet site consisting of the smallest side-chain residue glycine and hydrophobic residues (31IIGLMVGGVVI41). Binding energy analyses reveal the importance of π-stacking interactions, especially in the CHC site, hydrophobic interactions, and curvature matching. Strikingly, we find that the binding of DMF to the turn region can disrupt the D23-K28 salt-bridge that is important for the Aß fibrillation. These results provide molecular insight into the binding mechanism of fullerene to Aß protofibrils and offer new routes for the therapeutic drug design using fullerene derivatives against AD.

Transformation of an oblate-shaped vesicle induced by an adhering spherical particle.

We discuss vesicle-shaped transformations induced by an adhering spherical particle of different sizes using the classical Helfrich bending energy model. We focus on studying a model of an oblate vesicle with no spontaneous curvature in terms of fixed enclosing volume and surface area. Numerical solutions allow us to map out a complex phase diagram for such a simple, yet fundamental, physical system.